Exact zero-point energy shift in the e <g> (n E), t <g) (n H) many 
modes dynamic Jahn- Teller systems at strong coupling 

Nicola Manini* 
E.S.R.F., B.P. 220, F 38043 Grenoble Cedex, France 

Erio Tosattit 

International School for Advanced Studies (SISSA), via Beirut 2-4, 1-34014 Trieste, Italy 
Istituto Nazionale Fisica della Materia, Unitd Trieste SISS and 
International Centre for Theoretical Physics (ICTP), P.O. BOX 586, 1-34014 Trieste, Italy 

Abstract 

We find the exact semiclassical (strong coupling) zero-point energy shifts 
applicable to the e® {n E) and t<8> (n H) dynamic Jahn- Teller problems, for 
an arbitrary number n of discrete vibrational modes simultaneously coupled 
to one single electronic level. We also obtain an analytical formula for the 
frequency of the resulting normal modes, which has an attractive and ap- 
parently general Slater-Koster form. The limits of validity of this approach 
are assessed by comparison with O'Brien's previous effective- mode approach, 
and with accurate numerical diagonalizations. Numerical values obtained for 
t (8> (n H) with n = 8 and coupling constants appropriate to Cg are used for 
this purpose, and are discussed in the context of fullerene. 
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I. INTRODUCTION 



The dynamical J ami- Teller (DJT) effect, where one degenerate electronic level is nontriv- 
ially coupled to vibrational modes, represents perhaps the simplest type of problem where 
the Born-Oppenheimer approximation is invalid in principle, and a quantum mechanical 
treatment of the full electron-ion system is essential. Even for the simplest DJT case, true 
closed- form solutions have only been obtained for a restricted set of coupling values M. In 
spite of that, the DJT physics of a single vibrational mode coupled to the electronic state is 
quite well understood for all known symmetries, through analytic expansions valid in the 
opposite limits of strong and weak coupling, and otherwise through numerical diagonaliza- 
tion for intermediate values of the coupling. In a large number of cases of practical interest 
there are however many vibrational modes, coupled simultaneously to the same degenerate 
electronic state. This many- modes case has been extensively studied theoretically In 
particular, a first group of studies considers in detail the case of an impurity-related localized 
degenerate level interacting with the continuum of phonons in a crystal |||5| , while a second 
class of studies concentrates on the case of a finite, discrete set of vibrational modes [§]||[7|]. 
In the present paper we revisit the second type of problem, having in mind in particular the 
case of an isolated molecule/cluster /ion with a finite (although possibly large) number n of 
vibrational modes participating. 

To begin with, we recall the two special limits in which the many-modes problem is 
trivially solved 0. First, the case of all modes having the same symmetry and frequency, is 
simply reduced, by means of a rotation in the vibrational space, to the equivalent problem 
of a single mode coupled with the total coupling intensity, plus n — \ uncoupled modes. Sec- 
ondly, in the weak-coupling limit, degenerate perturbation theory applies, and contributions 
from different modes linearly superpose. Intermode interactions appear only at fourth and 
higher order in the couplings. 

The more realistic case of many modes with the same symmetry and different frequencies 
and couplings was apparently first addressed systematically by O'Brien [|J in the other 
customary limit, i.e., the strong-coupling semiclassic expansion. Her approach relies on the 
concept of replacing the n real modes with a single effective mode, perturbed by weak residual 
off-diagonal contributions. In this way one obtains an approximation which has been very 
useful, notably @ in the calculation of spectral shapes and reduction factors. Nonetheless, as 
we will show below, there is no well-defined limiting case where that approximation become 
exact in a controlled manner, except for the equal-frequency case. 

In this work we demonstrate an alternative approach, which is exact in the semiclassical 
(strong e-v coupling) limit, and applicable from strong to intermediate couplings. Rather 
than working in general, we address directly two specific and physically important cases 
of many- modes DJT, namely the linear e Cg> (nE) and t ® (n H). Our method relies on 



2 



a detailed analysis of the Born-Oppenheimer (BO) potential surface around its minimum, 
and of the normal frequencies of the classical small oscillations of the e-v coupled system. 
These frequencies, in turn, determine residual zero-point quantum correction to the classical 
motion on the BO surface. Comparison with the effective-mode approach M reveals that 
the latter, not taking into account these corrections, contains a systematic deviation in the 
strong-coupling limit, which we discuss in detail in Sec. |11 B| . Inclusion of these corrections 
generates the leading term of the exact strong-coupling expansion of the JT energy gain 
as a function of inverse coupling (Eq. (|16D ). As a very useful byproduct, we also obtain 
a detailed description of the low-energy vibronic spectrum, which is again accurate in the 
strong coupling limit. Finally, in order to provide a numerical application, we deliberately 
choose in Sec. |TTl|a borderline case of intermediate coupling, namely the negative ion of Ceo- 
That is a t <g> (n H) DJT problem, with n = 8, for which both exact and effective mode 
alternative calculations are feasible. Due to intermediate coupling, our approximation will 
not, of course, be particularly accurate. Nonetheless, being able to judge the sign and size 
of the deviation is especially useful in assessing errors, particularly in a case like this, which 
stretches a little beyond the reasonable borders of validity. 



II. THE E (8) (NE) SYSTEM 

We consider first the e <g> E linear JT Hamiltonian, a basic textbook p],|5]|| example of 
DJT, as well as a relevant model in many molecules and crystals. For convenience, we shall 
adopt the notations of Ref . Q . The Hamiltonian operator for the n-modes problem is 

h = \ E <* (p1 + «?) + E ( C L, 4) E f qa qi2 ) ( Cxa ) > (!) 

2 l V & ' i Vfts -qa) \CyaJ 

where u>i is the frequency and hi is the dimensionless coupling strength of the twofold- 
degenerate mode %. qi is a vector notation for the two normal coordinates qn and of mode i, 
and pi are the corresponding conjugate momenta, x and y label the two degenerate electronic 
states. The sums should be understood as J27=i- Note that we use the second-quantized 
notation for the fermions, and the coordinate description for the vibrational degrees of 
freedom. We also set h — 1, thus making no distinctions between (angular) frequencies and 
energies. 

As a first step, we treat the q} as classical coordinates, and study the (lowest) Born- 
Oppenheimer potential surface: 

v (a) > = E^(ft) 2 + min E kiLJi (W4 c z - c l c v) + ft2(4 c s/ + c l c x)) W) (2) 

i % 

The minimum over the fermionic degree of freedom is the lowest eigenvalue of 
the 2x2 electronic problem representing one fermion in the degenerate level, which is 
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1 /2 

— (j^i ki^iQi + 2 J2i<j hkjUJiUJjqi ■ qA . The lowest BO surface is conveniently rewritten 
in polar coordinates q\ = cos 9 iy qi sin 0j) as 



1/2 

^ (?i> 0<) = g II ^ ~ I ? fe < w * 2 5? + 2 £ kikjUJiUjqiqj cos(^ - [• (3) 



In this form, the minima of V, corresponding to the classical stable equilibrium configura- 
tions of the system, are straightforwardly discussed. In order to minimize V, the argument 
of the square root should be maximum, which is obtained when all 9i — 8j simultaneously 
vanish. In this case, for all cosines equal to one, the square root can be explicitly executed, to 
obtain separate dependences on the different q^s. The minimum of V is therefore obtained 
for 

^min _ t 
Hi — K i 

ef n = e, (4) 

where 9 is arbitrary. As form (§) explicitly shows, the potential energy is independent of 
the common rotation angle 9 = J2i@i/ n - This implies in particular that the configurations 
of minimum potential energy constitute a continuous manifold, parametrized by 9, with the 
topology of a circle, exactly the same as in the one-mode problem. This result was to be 
expected also on the basis of more abstract considerations 0, and we shall return to it 
below, when analyzing the small oscillations around the minimum. 

The value of the potential at the minimum is the classical JT energy gain 

^class = -~ k i U i = -^effWeff , (5) 

where we have introduced, following Ref. ||, = J2i kf and u e fi = J2i kfuii/kl s . 

This was just a re-derivation of well-known results. The next step is to compute the 
quantum corrections to this classical result. The quantum corrections add to -E c iass in the 
form of a strong-coupling expansion in negative powers of k e s. The leading term, of order 
zero, is a shift independent of k eS , due to the change of zero-point energy of the system 
following the JT coupling. 

This is most easily illustrated in the case of just one mode. The ground-state energy 
of Eq. ([I]), E{k\ = 0) (no JT coupling), is just the zero-point energy of the (two-fold 
degenerate) mode, amounting to 2 • On the other hand, when the coupling is very 

strong, the dynamics factorizes into a radial harmonic mode plus a free pseudo-rotation. 
The harmonic zero-point energy is now only For a single mode, the total JT energy 

gain is therefore, as is well known, 

E JT - E(ki = 0) = E class - iwi + -wifcrV + { k i 4 ) > ( 6 ) 
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where -E c i ass = — \k\uii term represents the lowering of the BO potential minimum, —\uJ\ 
representing the zero-point energy gain, and the j(= ±~ for the ground state) term is the 
residual zero-point energy associated with quantization of the (pseudo-rotational) 9 motion. 
In the rest of this section we wish to generalize the second and third term in Eq. to the 
many-modes case. 



A. The semiclassical expansion 



We now extend this kind of semiclassical expansion to the many-modes case. It is clear 
that the zero-point energy for zero coupling E(ki — 0) = 2 x \ To determine the 

zero-point energy for the motion around the many-modes potential minimum in the finite 
coupling case, we expand the BO potential to second order about the minimum in the 
x = (qi, q 2 , 9i, 8 2 , ■■■) coordinates, and compute the classical normal modes of vibration: 



^ 2n 

V (x) = E dass + - 

A* if 



d 2 V 



x, 



doc n dec v 



x: 



lin ) + O (( 
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where the Hessian matrix of the derivatives is 



V 



and E = (£i kf uji) . 

The kinetic energy for these coordinates is 
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The normal frequencies uo 2 and vibronic modes are obtained as eigenvalues and eigenvectors 
of the matrix 



T (x min j • V 



V 



-I 











ujf — Eklujf —Ek 2 u 2 uf ... 





—EkfuiU% uj\ — Ek\u)\ ... 



(11) 



This matrix is non Hermitian, but has real eigenvalues. The block-diagonal form shows that 
the small oscillations of the radial and angular variables are uncoupled. Each radial coordi- 
nate qi corresponds directly to one radial mode of the same frequency Ui as the uncoupled 
modes. On the contrary, the dynamics of the angular variables are intercoupled through 
the off-diagonal elements in Eq. flTTl). The secular equation for the angular eigenvector of 
components Xi and eigenfrequency u 2 is 



V 2 Xi = ojfxi 



jUjXj 



(12) 



with the solution 



X;, = const 



^'7 



UJ? — u 2 



(13) 



where the corresponding v 2 is a solution of the equation 

k 2 t,) 3 



(14) 



As a simple graphical analysis suggests, this equation has as many solutions v 2 as original 
modes. For any frequencies Uj and couplings kj, the lowest solution is V\ — 0, corresponding 
to an eigenvector Xi = (1,1,1,...), i.e. to the totally symmetrical coordinate {J2@i/ n )'- it 
represents the "soft mode" of pseudo-rotation along the circular JT valley. If the other 
solutions uj are sorted in ascending order of frequency, the inequality Uj-\ < Uj < ujj can 
be seen to hold. In the special case of modes of the same frequency u)j-\ = ujj, the implied 
equality Uj = ujj is also true. Indeed, a set of d modes with equal frequencies can be 
rewritten as a single coupled mode, plus d — 1 uncoupled ones: the unchanged frequencies 
Uj = Uj correspond therefore to the angular part of the uncoupled modes. In the simple, 
but instructive, case of n = 2 modes, the (single) nonzero angular frequency is 

1/2 



U 2 



U\<j)2{k 2 u)\ + k x u) 2 ) I {k^uJx + k 2 uj 2 ) 



(15) 



6 



This expression for the new angular frequency has the form of a weighted geometrical mean 
between the unperturbed frequencies, with such weights that v\ is attracted towards the 
frequency of the mode with weaker coupling. This is reasonable because the frequency of 
any mode whose coupling is exactly zero should of course remain exactly unchanged. The 
frequencies Uj display the same behavior in the general case: the new modes Vj are located 
in the intervals between neighboring frequencies uij-i, uuj, attracted towards the mode with 
weaker coupling strength kj. In the limiting case kj = 0, therefore, Vj restores smoothly the 
twofold-degeneracy of the uncoupled mode Uj. For a point impurity coupled to the phonon 
continuum, equations analogous to Eq. ([14]) were derived earlier || , but apparently not put 
to practical use. Here, we will use our eigenvalue equation ([14]) for the calculation of vibron 
frequencies, and of the associated zero-point energy. 

We underline, incidentally, that this equation has a very familiar form, namely that of a 
Slater-Koster scattering problem for a separable potential, or of a collective-mode equation 
|T0| , or of a BCS gap equation [pp . By coupling to the same electronic state, the different 
vibrational modes get effectively coupled to each other by a kind of separable "attractive 
5-function" potential. The lowest bound state is forced to zero frequency, v\ = 0, by 
what appears to be an exact sum rule, related to the 0(2) symmetry of the Hamiltonian 
(0), corresponding to the independence of the energy of the angle 6. In the BCS case, 
broken gauge invariance and Goldstone's theorem give rise to a formally similar sum rule. A 
universal feature of linear JT systems, the zero mode is also present in the phonon continuum 
case 0- 

In the evaluation of the zero-point energy, the normal angular frequencies appear only 
as a sum: 

£ JT - E{h = o) = + M£"i + E fi ) - 1> + °(Ki) 

V i angular modes Uj J i 

= ~\k 2 cS uj cS + 1 - Uj ) + 0(Ki) . (16) 

This main result is not in "closed form" , since the uj are defined implicitly as the solutions 
of Eq. (|Hp . However this is a very marginal shortcoming, involving in each specific case 
the simple solution of a numerical equation, explicitly given in terms of the parameters uij 
and kj. We also remark that a rescaling of all the couplings kj — > akj brings no change to 
the secular equation flT2"p: the angular frequencies Uj, therefore, do not depend on k e g, but 
only on the ratios among the couplings. Thus, in particular, the v j ~ zero-point term 
is really a term of order zero in k e g, as it should. 

This approach provides quite naturally the detailed structure of the spectrum of the ex- 
cited states in the strong-coupling many-modes case. All the uncoupled frequencies Ui persist 
as radial frequencies in the DJT spectrum, while new harmonic angular modes of interme- 
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diate frequencies Uj (solutions of Eq. (|14D) appear in between, as illustrated schematically 
in Fig. 0. 

In addition to these main spectral structure, the zero-frequency mode v\ introduces an 
entire ladder of low-energy excitations, corresponding to the quantization of the pseudo- 
rotation around the manifold of minimum BO potential. We write these excitations in 
terms of the "moment of inertia" of the pseudo-rotor associated with the J2 ®% coordinate 

i- 

e "W^-»' (17) 

where the allowed values of m are m = ±|, ±|, as required by the Berry-phase prescrip- 
tion. E rot (m) adds to Ejt to give the low-energy states. In particular, the energy of the 
(twofold degenerate) ground state is given in this scheme by 



? GS « £jt + E rot ( ±^ ) . (18) 



B. Numerical calculations for e (g> (nE) 

The relative accuracy of the results just derived and of O'Brien's effective-mode approx- 
imation m 

Efl s - EMh = 0) = -kVcff - k ff - ^ 2 ~^ 2 + O(kJ) 

i i 4 cu cS 

= -i^ s -i^i^ + Ofc 2 ) (19) 

L 4 LU cS 

where uj 2 = J2i kfujf/kl^, can be assessed by comparison with the exact numerical diago- 
nalization results of ([I]) on a truncated oscillator basis. The ground-state energy obtained 
by numerical diagonalization is a variational estimate of the exact ground-state energy. It 
converges very rapidly as the number N of included oscillator states exceeds k 2 s . In practice, 
for k 2 s ~ 100, inclusion of iV = 100 quanta in the oscillators ladder yields an accuracy of 
order 10~ 8 u;cfr, largely sufficient for our purposes. 

For the case of only two modes, where an explicit form (|15D is available for i>2, the 
difference between the two is 

E?l s - E JT = - ^ 2 + ^ cff2 + + uj 2 - u 2 ) + O(kJ) . (20) 

4 LU c g I 

As a first example, we consider the case of two modes with similar frequencies, such that 
co>2 — ^i = Su is a small parameter. The difference in ground-state energy give by the two 
approximations is 



S 



^" - E " « f + 2(kl - k\)(k^f ( -J^V + ... (21) 

This difference is therefore very small for small 5u: in this limit the two expressions are 
essentially coincident, and we verified that they both agree with the exact numerical ground- 
state energy at strong coupling. This is not surprising, since we know that the effective- 
mode result is exact in the trivial case of equal frequencies ||, and it confirms that also our 
expansion is correct in this limit. 

Moving to the more interesting situation of very different frequencies, for example uj 2 = 
10 • Ui, we find that the correction becomes more important (~ 5% oo c g). Fig. ||] shows the 
error of the approximate expressions AE = _E a PP rox — ^exact f Qr ^ e ground-state energy 
at different values of the coupling. Our approximate formula flT6|) converges systematically 
(from below) to the exact energy modulo corrections ~ k~$ . The effective mode expression 
(|IT5|) differs from Eq. (^) for a quantity depending on the individual frequencies and on 
the ratio between the couplings, but not on the total coupling strength k~$. This difference 
introduces a systematic shift in the strong-coupling limit of the the expression ( fL9|) . 

r -1 -2 

In the extreme case of very small frequency ratio U1/U2 (<< f~ + )> the systematic 
shift introduced by Eq. fllTf) becomes as large as 



which is a relevant fraction of u e g. 





~ 2 ) + 2 




^ + .-, (22) 
uj 2 



A third example (Fig. |3|) illustrates the validity of the method for the case of three 
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modes. For a rather large frequency spread - ~^ cff (~ 2), the approximate formula (16) still 
converges to the exact energy, but not monotonically in k e g. 

An interesting observation suggested by Figs. ^| and |3] is that, for weak enough coupling 
- say k cS < 3, the effective-mode theory, owing to the systematic shift discussed above, 
can yield better agreement with the exact ground-state energy than our method, which is 
instead superior at strong coupling. Thus the effective-mode method and the present one 
are to some extent complementary. 

Besides the ground state, the Lanczos technique employed for diagonalization easily 
generates a few low-lying excited states. These accurate excitation energies can be compared 
with our approximate frequencies UiS. In Fig. ^ we report the excitation energy of the low- 
lying m=±| states. These energies show a clear convergence to the fundamentals, overtones 
and combination states of uj\ and as expected from the theory. Above each of these 
origins, higher-m states form pseudorotational ladders, as determined by E mt (m). In Fig. |] 
we represent for clarity only the m = ±| excitation above the ground state, and compare 
it with the theoretical value given by Eq. (|17|) . Interestingly, a parallel comparison with the 
effective-mode model shows that, even though the ground-state energy estimates in the two 
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models differ, the values of the pseudorotational quantum of energy (expressed in completely 
different forms) are numerically coincident ||12|| . 

We stress that our adiabatic potential (fj) does not explicitly include centrifugal terms, 
unlike what is customarily done away from strong coupling. Reference 0], for example, de- 
scribes a self-consistent prescription for embodying these terms at the outset. For the sake of 
completeness we include in Fig. |] the excitation energy to the first m — ±| pseudorotational 
state obtained through that method. The self-consistent procedure renormalizes the inertial 
moment, eliminating the k~$ divergence at weak coupling, but it ends up providing a worse 
approximation to the exact energy in the region of intermediate to strong coupling. Hence, 
inclusion of centrifugal terms does not appear to be useful in this regime. Concerning the 
0(k~ff) terms, we note that, while in Fig. |^ the pseudorotational contribution E Iot f±§) 
shifts the approximate ground-state energy closer to the exact one, in Fig. | for k e s > 4 the 
error of the approximate expression (|T6|) is positive: thus addition of the pseudorotational 
contribution moves the approximation further away from the correct value. We conclude 
that the pseudorotational contribution as given by Eq. (|TTD, contrary to the single-mode 
case, does not exhaust all the 0(k~^) corrections to the truncated expansion (|T6|). We will 
not further investigate such terms in the present work. 



III. MANY MODES T ® (N H) 



After the basic e®(nE) system, treated in the previous sections, we wish now to consider 
a second, different case, to illustrate how much of the procedure used can be carried over. 
The interaction of an orbital triplet (t) with a set of fivefold vibrational modes (H icosahedral 
representations) is our next choice. As we shall see, apart from some technical differences, 
we will be able to follow very closely the approach that proved successful for the e ® (n E) 
case. 



We start with the Hamiltonian operator fl3| . [T4 
1 



H 



(23) 





( Qa - V3qi4 


-V3q i3 




A 


( c xa \ 


C \ai C ycn C za) ^2 


-V3q i3 


q a + \/3gi4 






C y(7 


i 


{ -V3q i2 






j 


\c za ) 



where now the three degenerate electronic states are labeled x, y, z, and ^ indicate a five- 
dimensional vector of components (qn, ...qis). Reference |14| introduces a polar parametriza- 
tion of the five-dimensional space of one single mode i in term of a radial coordinate qi, plus 
four angles c^, 9i, (pi, and ji. The same work |TJ| reports for a single mode i, the expression 



for the lowest electronic eigenvalue of the one-electron matrix in terms of these coordinates: 
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- cos(ai)kiUiqi 



(24) 



The BO potential (obtained adding to the electronic part) is therefore independent 

of 9i, 4>i, and but it depends explicitly on the angle aj. In particular, it is minimum at 
qi = ki, da = 0. For this special value, the angular parametrization of q} becomes singular, 
the coordinate q\ being independent of the angle 7$. The DJT valley of t (g) (n H) is thus, as 
well known, two-dimensional, parametrized by 9i, (pi. 

Now, let us move on to n > 1 fivefold modes. Consider the vector 5Z i kiUJiqi intervening 
in the JT coupling matrix (|23|) : if we indicate with a, 9, (ft, and 7 the corresponding set of 
polar angles, the lowest electronic eigenvalue for the general case writes 



cos m 



k i u S 



1/2 

cos(a<) ( ^i^ili + 2 kikjUJiUJjqi ■ q$ ] . (25) 

i<j 



The minimization of the corresponding the BO potential, function of 5n variables, is therefore 
very similar to the analogous for Eq. @: here the electronic term is expressed as the negative 
of a product, where both terms can be maximized at the same time. As a result, at the BO 
minimum, a = 0, and all the q^ should align with one another in a common direction: 

^min _ 7 
q% — K i 

af n = a = 

6f n = 9 (26) 
C in = 

7f n = 7 

Again the trough is two dimensional, parametrized by 9, <ft, since 7 is singular as discussed 
in the one mode case above. The expression (H) for the classical JT energy gain holds in the 
t ® (n H) case. 

The next step is to generate the quantum corrections to the classical result. For this 
purpose we need the harmonic frequencies of oscillation around the minimum (|26|) . Here a 
completely analytical approach, as in Sec. |H, fails, for two reasons. As a first point, the factor 
cos(a), contributing to the Hessian matrix 9 ® 2 q x , is unavailable as an explicit expression in 
terms of the coordinates x = (qi, q 2 , a^, a 2 , 0i, ■■■■)■ Secondly, the parametrization is 
singular right at the minimum, thus it cannot generate all the angular modes. 

We resort therefore to an alternative, more numerical approach. We express explicitly 
the lowest electronic eigenvalue in Cartesian coordinates x = (qn, q^, q 2 i, fe) com- 
pute the 5n x 5n Hessian matrix (at an arbitrary minimum point in the trough) and the 
kinetic matrix [see Eq. (|?p, (Q)], and diagonalize the product T (a? 11111 ) • V [see Eq. (PI)]. 
For n — 2, the calculation can be done analytically, choosing some special points in the 
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minimum manifold (p6|) , such as 6 = 0, <fi =anything [i.e. <fj = (^-,0,0,0,0)], or 9 = tt/4, 
= 7r/2 [i.e. gi = (fcj/4, 0, 0, —^1/3/4, fci\/3/2)]. The resulting matrix has two threefold- 
degenerate eigenvalues of frequency u% and u 2 respectively, a twofold eigenvalue v x = 0, and 
a second twofold eigenvalue of frequency z/ 2 given by Eq. (|15|). The eigenvalues are of course 
independent of the choice of the minimum point around which the expansion is done. In the 
general case of larger n, we diagonalized numerically the dynamical matrix, obtaining for 
any choice of the frequencies and couplings n threefold-degenerate eigenvalues ui,... u n , and 
n new twofold-degenerate eigenvalues, which we called Vi. Again we have always v\ = 0, and 
all z/j's are located between subsequent modes u^i and u^. Thus, directed by the similarity 
with the e®{n E) case, we verified by substitution that the new frequencies z/j are solutions 
of Eq. (0). 

From the above analysis we conclude therefore that (i) in analogy to the e®(nE) case, 
each original frequency u>i is still present in the e-v coupled spectrum; (ii) in t <g> (n H) 
each uji corresponds not just to the radial degrees of freedom, but includes pairs of coupled 
angular modes as well, for a total degeneracy of three; (iii) new harmonic vibron modes 
appear at the same frequencies - dictated by the same equation ([TJ|) - as in the e ® (n E) 
case (see Fig. [I]); (iv) in the t®{nH) case, the new modes z/j are all twofold degenerate; (v) 
in particular, the lowest new frequency ui = again corresponds to the free pseudorotation 
around the minimum trough; (vi) finally, the semiclassical DJT energy gain, in analogy to 
Eq. (|16|), is for this case: 

E 3T - E{ki = 0) = -^ c Voff + ~ I 3 ]T ^ + 2 J2 I - \ E "i + °(Ki) 

= -^eVeff + $>i - Uj) + 0(k;i) . (27) 

i 



A. C 60 : a case of intermediate coupling t <g> (n H) 

As an application, chosen to test the limits of applicability of our approximation, we 
consider the DJT problem of the fullerene anion |IJ],[TJ|-|nj] , where a single electron in a ti u 
electronic state couples to n = 8 H g vibrational modes. The frequencies and couplings of the 
eight modes differ strongly, and none of the couplings is particularly large. The numerical 
values of the parameters for the Cg ion, the same adopted in our previous work |L7|, are 



reported in Table |. The table gives also the new vibron frequencies uj of the the coupled 
system, calculated within our approximation. Note that mode 6 does not give rise to a 
new frequency, since its coupling is vanishingly small. On the basis on these numbers, we 
have computed the ground-state energy gain in the approximation (p7|), and we report it in 
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Table H, along with the corresponding values in the effective-mode scheme, as presented in 
Ref. JUj, and with the exact result, obtained by Lanczos diagonalization [|17| .|T8|1 . 

We note that the two approximations, the effective mode and ours, are essentially equiv- 
alent, both of them in error by some 25% relative to the exact result. The approximate 
energy gains, in particular, correct in excess the initial classical JT gain |_E c iass| (too small, 
by about a factor 1/2) and are now 25% too large. The reason for this is that the vibron 
zero-point energy, here larger than the classical energy gain itself, overestimates the true 
quantum correction, indicating important higher-order corrections. This reflects the fact 
that, in Cg , & e fr is only slightly larger than unity, since the individual couplings are rather 
weak, so that a semiclassical expansion truncated omitting terms of 0(/c e ^ r 2 )(~ 250cm~ 1 ) and 
higher is quite far from convergence. Qualitatively, this regime corresponds to the region 
k 2 eS « 2.5 in Fig. |. 

In the detail, the effective-mode energy is a few wave numbers better than our approxi- 
mation. The reason for that is the asymptotic shift previously discussed for the e <S> (n E) 
case, and illustrated in Figs. [] and ^ which makes the effective mode competitive when the 
coupling is weaker. 

If in addition we consider the contribution, of order k~$, of pseudo- rotations, we can 
translate Eq. ([L7]) to: 

E™\L) = —^- I .L(L + l). (28) 
6 £ kfuJi 



The inclusion of this contribution (60. lcm -1 ) for the L = 1 ground state jl3],|14] brings the 
semiclassical estimate of the ground-state energy gain to within 20% of the exact numerical 
value (Table [H]). Note however that, according to (|28"D, the first-excited (L = 3) pseudorota- 
tional state should be found at ~300cm _1 above the ground state, that is slightly above u)\. 
This confirms that the fullerene ion is really in an intermediate-coupling situation, rather 
far from all approximate limiting regimes. 



IV. DISCUSSION 



O'Brien's approach - based on an "effective- mode" picture - is available in the literature 
for the description of the low-energy states of a many-modes JT system. It introduces 
a fictitious, effective single mode, treating perturbatively the residual corrections. This 
approach, conceptually attractive as it is (and exact in the limit of equal frequencies Ui = u) 
H, is not completely satisfactory, particularly when the spread of the frequencies is sizable, 
and couplings are large. In that case, the effective-mode expression (|19]) for the ground-state 
energy gain has a deviation which may be a large fraction of the effective frequency uj c q. This 
shift, for given frequencies and ratios among the couplings, is independent of k e g, therefore 
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the inclusion of higher power corrections in k~$ to Eq. ([19]), such as those of Eq. (44) of 
Ref. || cannot correct the asymptotic behavior at large k e g. 

We have introduced a semiclassical treatment of the many-modes problem, which corrects 
the difficulties of the effective-mode model. We have applied this method to two specific JT 
cases making it clear that, with straightforward modifications, it can be extended to other 
D JT systems. The method is based on the approximate calculation of the new vibron normal 
mode frequencies v\ arising after DJT coupling. Of these modes, n — 1 fall in the intervals 
between successive original bare mode frequencies, the remaining one is a zero mode V\ = 0. 
The approximate z/j's are again solutions of a collective-mode-type equation [fLOfl , already 
obtained in an equivalent form in previous studies of the continuum case ||. Numerical 
tests confirm this picture, including the recovery of the original, unperturbed frequencies uj^s 
in the strong-coupling spectrum. 

In particular, we obtain expressions ([T6|j27|) , in terms of these frequencies, for the ground- 
state energy which are exact in the strong-coupling limit. At finite coupling, our formulas 
introduce an error of the order k~g, which is only partly corrected by the term Eqs. (|T7|) , (|28|) 
of pseudorotation along the BO minimum trough with the associated Berry phase constraints 
|13| . However, in cases with strong enough coupling for our expansion to hold for the 
potential, the expressions to order 0(k® s ) gives satisfactory results, and will not need further 
refinement. 

We provide ample numerical verification for our expectations in the e ® (n E) case. A 
detailed numerical test of the result for the case of t®(n H) would be rather difficult, because 
of the explosion of the basis set size, which occurs when many vibron states are included, as 
necessary at strong coupling. The case of Cg , not a very strongly-coupled system, therefore 
amenable to exact diagonalization, shows that our approach provides results for the energy 
gain which are basically equivalent to those of the effective-mode approximation. 

More interesting would be a direct experimental confirmation of the new vibron frequen- 
cies which we have found. The ideal test system could be a molecular system where the 
couplings to a few (2 or 3) modes are known, reasonably strong, controlled, and possibly 
switchable on and off (for example upon doping): in such a situation, infrared/Raman spec- 
troscopy should easily evidence the birth of new modes. Benzenoid cations might represent 
a suitable systems for such future investigations. 

The present formulation of course neglects a large set of effects potentially capable of 
modifying the picture to some degree, especially on the low- and high-energy ends of the 
vibronic spectra. The most obvious limitation regards the assumed harmonicity of the "small 
oscillations" (^), and the neglect from the beginning of nonlinear terms in ([I]). In fact, there 
is no truly flat JT valley in a real system. The effect of valley warping is to localize the 
pseudorotational motion into a discrete set of minima. The importance of this phenomenon 
is determined by the relative value of the pseudorotational energy E rot and the height of 
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the barriers between minima. Initially the warping will affect mostly the zero-energy mode, 
but for very large JT distortions, the whole shape of the BO surface will change. As a 
consequence, our results apply best to a strong-coupling case where the large value of k e g is 
given by the addition of many individually small contributions fcj, such that the distortion 
of each mode's coordinate remains small, and the effect of higher-order terms is therefore 
weak. 
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TABLES 



i 


Wi(cm _1 ) 


ki 


^kfuji(cm x ) 


z^(cm _1 ) 


1 


270.0 


0.868 


101.71 


0.0 


2 


430.5 


0.924 


183.78 


329.5 


3 


708.5 


0.405 


58.11 


633.8 


4 


772.5 


0.448 


77.52 


742.3 


5 


1099.0 


0.325 


58.04 


1031.8 


6 


1248.0 


0.000 


0.00 


1248.0 


7 


1426.0 


0.368 


96.56 


1302.4 


8 


1575.0 


0.368 


106.65 


1519.0 




(Jeff = 581.1 


k cS = 1.532 


£ciass = E s*f = 682.36 





TABLE I. The frequency /coupling parameters for many-modes DJT of C^ . The last line 
combines them to give the effective parameters introduced in Ref. With k^ s = 2.35, Cg 
is indeed an intermediate-coupling system. The last column reports the new frequencies of the 
DJT-coupled system. 



DJT energy gain 



this work 



Ref. HI! | 



exact [17] 



order k^ s 
order k® s 



pseudorotation k eS term 
Total gain (no pseudorotation) 
Total gain (including pseudorotation) 



\ E kfuJi 



-682.4 



the same 



; ujj) = -722.8 



^r 1 



-i 



60.1 



-1405.1 
-1345.0 



2 u eS 

the same 
-1395.8 
-1335.7 



713.4 



-1125.7 



TABLE II. The relevant ground-state energetics (in cm -1 ) for the many-modes DJT of Cg , in 
the present model and as obtained in the cited works by O'Brien. The classical potential-lowering 
term, the lowest order quantum corrections, and the pseudorotational contribution are listed. 
Finally, the strong-coupling energy gains (up to order k® s , and including the pseudorotational 
correction) are compared to the result of Lanczos diagonalization. 
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FIGURES 



OLD FREQUENCIES NEW 




FIG. 1. A schematic picture of the spectrum of vibronic levels of a coupled many-modes JT 
system, in the strong-coupling limit. For simplicity, only the harmonic one-phonon (fundamental) 
states are drawn. The thin lines represent a few low-lying pseudorotational levels. In the e®{nE) 
case the harmonic vibrations are nondegenerate, but all the states acquire a twofold degeneracy 
when the pseudorotation is considered. The u)i and V{ fundamentals in t <S> (n H) are threefold and 
twofold degenerate, respectively, and these degeneracies combine further with those of the "soft" 
states of pseudorotation around the trough. 
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30 60 

k 2 

^eff 

FIG. 2. The difference between the approximate expressions and the exact ground-state 
energy, obtained by Lanczos diagonalization including up to N = 100 oscillator states, for two 
modes of frequencies ui2 = 10 • ui\ , as a function of the total JT coupling strength k^ = 2k\ = 2k\. 
The figure shows the semiclassical expression its correction including the pseudo-rotation 

contribution (|l^), and the effective- mode formula (|I^), with also the complete version - Eq. (44) 
of Ref. - including the k~^ correction. 
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0.2 



3° 0.0 



-0.2 



-0.4 








10 



effective mode Eq.(19) 
effective mode to 0(k ef{ ) 



present work Eq.(16,18) 



20 



30 



40 



v eff 



FIG. 3. The difference between the approximate expressions and the exact ground-state 
energy, obtained by Lanczos diagonalization including up to N = 40 oscillator states, for n = 3 
modes of frequencies uj\ = W2/IOO = U3/I2O, and couplings k\ = 2/C2 = 2^3, as a function of k^ s . 



The semiclassical expression (16) and its correction including the pseudo-rotation contribution 
(|l8|) are indistinguishable on the scale of the figure. The figure shows also AE for the perturbative 



formula (|19|), and for the version including the k eS correction. 
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3 




1 1 1 1 1 1 u 

10 20 30 




FIG. 4. Excitation energy (solid lines) of selected low-lying excited states for n = 2 modes of 
frequencies oj\ = 1, 002 = 3, and coupling ratio k\ = 5^2 (solid lines). The four lowest m = ±5 
states and the lowest m = ±| state are drawn, and compared with the theoretical previsions 
(dashed lines) and assignments of our model, valid at strong coupling. The dotted line represents 
the excitation energy of the lowest m = ±| state including the self-consistent centrifugal correction 
of Ref. @. 
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